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I. INTRODUCTION 



Chemical reactions of the early actinide elements in aqueous solution are complex and challenging to predict. 
Elements U, Np, and Pu each have four or more oxidation states in acidic environments (Am has three). Two of these 
states (III and IV) are hydrated An 3+ and An 4+ ions; the other two (V and VI) form linear trans-dioxo AnOj and 



An0 2+ actinyl complexes. Disproportionation reactions are common, especially in the case of plutonium for which 
three of the redox potentials are nearly the same (about one volt). As different oxidation states have widely different 
solubilities 1 ', redox chemistry plays a crucial role in the environmental dispersal of actinidesP The rich behavior 
of actinide ions may be traced to the valence electrons, especially those in the 5f shellP Because the 5f orbitals 
are relatively localized, the Coulomb interaction induces strong correlations between electrons in the shell. The 
possibility that strong correlations in these solvated actinide and actinyl ions enhance tendencies to disproportionate 
is an intriguing hypothesis. Support for this idea is provided by a Hubbard model of 5f electrons that disproportionates 
when solved in the Hartree-Fock (HF) approximations' 

Quantitative models of actinide reactions must overcome several obstaclesP First relativity and the energetics 
of solvation must be taken into account. At present only the density functional theory (DFT) method is capable 
of modeling these aspects accurately. However, incorporating the physics of strong electronic correlations among 
the 5f electrons presents a greater challenge as these are known to be poorly captured by DFTP In this paper 
we take a hybrid approach to solving these problems by using DFT to construct a generalized Anderson impurity 
model of the frontier orbitalsP^^ Exact diagonalization of the Anderson impurity model corrects the free energy 
obtained from DFT alone, yielding improved predictions for redox free energies. We emphasize that the hybrid 
method outlined here differs markedly from the LDA+U approach as it does not simply modify the LDA functional 
to partly account for the Coulomb repulsion; rather a high-dimensional many-electron Hamiltonian that models the 
physics of strong correlations between the important low-energy states is diagonalized. The approach also differs 
from configuration-interaction (CI) and its variants^ in two significant ways: First we are able to exactly diagonalize 
the low-energy effective model with no restrictions placed on the ground state wavefunction. Second, the two-body 
Coulomb interaction is not the bare electron-electron repulsion but rather an effective interaction that takes into 
account the effects of screening. 

The rest of the paper proceeds as follows: In Sec. [H]we discuss the ab initio part of the calculation and compare 
the results we obtain to calculations by other workers. The independent-particle model that describes the Kohn- 
Sham (KS) orbitals and the spin-orbit interaction is introduced in Sec. |III| Incorporation of the effective Coulomb 



interaction via a many-body model of the low-energy degrees of freedom is carried out in Sec. IV The physics of 
strong correlations are illustrated with the use of a simplified model in Sec. [V] The calculation of redox potentials 
and other observables by the hybrid approach are presented in Sec. |VI| We conclude with some discussion in Sec. 

eh 



II. DENSITY FUNCTIONAL THEORY 



DFT based studies of early actinides in aqueous solution have modeled the structure, vibrational frequencies, and 
free energies of hydrationP3H251We j n0 ^ attempt to review this work comprehensively here but refer the reader instead 
to Ref. 5 and references therein. As there is significant hybridization between the higher orbitals, accurate calcul ations 
require full quantum mechanical treatment of all orbitals with principal quantum number n = 5 and higher! 26 ' I Using 
small cores, ab initio calculations performed by Shamov and SchreckenbacfP^ were able to reproduce An(VI)/An(V) 
redox potentials to within 0.6 volt of the measured values. 

The Amsterdam Density Functional (ADF)PS1 is an attractive DFT package for modeling actinides in solution 
because it includes relativistic corrections via the zeroth-order regular approximation (ZORA)- 29 ' 30 !, uses a basis of 
localized Slater-type atomic orbitals, and models solvation with the Conductor like Screening Model (COSMO) In 
our calculations the first coordination sphere of water molecules are treated quantum mechanically; COSMO is used to 
simulate a bulk dielectric medium beyond the sphere. Quantum mechanical modeling of the second sphere of hydration 
may lead to improved agreement with experiment)^ but we defer that for future work. For the sake of simplicity the 
number of water molecules in the first solvation sphere are kept constant for each oxidation state: Eight each in the 
case of An(III) and (IV) and five for An02 (V) and (VI) . (As discussed later, we find that coordinating Pu(III) with 9 
water molecules makes only a small difference.) We use the revised PBE exchange-correlation functional E2HHp or the 
actinides we employ a basis set modified from ADF's triple-C, doubly polarized (TZ2P) ZORA wavef unctions. The 
frozen core consists of the 60 orbitals for which n < 4. The ADF-supplied basis set has 78 frozen core orbitals; 60 of 
these are retained in the core and the 5s, 5p, and 5d orbitals are promoted to valence orbitals. The number of sets of 
fit functions is accordingly increased from 79 to 87. For oxygen and hydrogen the relativistic TZ2P all-electron basis 
set provided by ADF is employed with no modification. 
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(a) Pu(H 2 0)!+(m) (b) Pu(H 2 0)^+(IV) 





(c) Pu0 2 (H 2 0)+(V) (d) Pu0 2 (H 2 0)2+(VI) 

FIG. 1. Converged geometries of the plutonium complexes. 

Spin-unrestricted calculations are carried out in three stages. As a first step the geometry of the actinide (or actinyl) 
plus the first solvation sphere is optimized in the gas-phase (no COSMO). We verify that all vibrational modes have 
real frequencies at the optimal geometry; these frequencies are used later in the thermodynamic calculations. The 
geometry is then allowed to relax in solvation within the COSMO approximation for the surrounding dielectric, with 
the following cavity radii: 1.350 A (H), 1.517 A (O) and 2.10 A (An). Thermodynamic properties are calculated 
in the ideal gas approximation with an effective pressure of 1354 atmospheres to account for the reduced entropy of 
translation as appropriate for the aqueous environment PSl As a final step an independent-particle model of the frontier 



orbitals is constructed (see Sec. Ill) and the single-electron spin-orbit interaction is added to correct the electronic 
contribution to the free energy. 

Minimum energy configurations are found to have highly symmetrical geometries consistent with X-ray absorption 
fine structure (XAFS) and extended- XAFS (EXAFS) measurements. For example, the geometry of the Pu(III) and 
Pu(IV) complexes converge to the cubic configuration shown in Fig. nlwliile the actinyls form pentagonal bipyramids. 
Interestingly, the converged geometry of the Pu(V) actinyl complex (Fig. [I] (c)) rotates one of the water molecules so 
that both hydrogen atoms are on the same side of the molecular equatorial plane. This distortion of the molecular 
symmetry is also seen in the converg ed geometries of U(VI), Np(VI), Pu(V), Am(V) and Am(VI). Table [I] shows 
the calculated bond len gths of the four oxidation states. The calculated geometries compare well with those found 
in previous DFT studies ^ * 18 * 21 * ^ of actinyls in solution and also, as shown in Fig. [2] with available experimental 
measurements, though DFT overestimates the mean An-OH2 bond lengths by up to 0.1 A in the case of the actinyls. 

To calculate the redox potentials we follow the procedure outlined in Refs. [17] and [24J, extending that work to 
include the An(V)/An(IV) and An(IV)/An(III) redox couples as well as reactions involving americium. The three 
half-reactions are as follows: 

An0 2 (H 2 0) 2+ (VI) +e- -> An0 2 (H 2 0)+(V) 
An0 2 (H 2 0)+(V) + 4H 3 0+ + e" -> An(H 2 0)*+(IV) + 3H 2 

An(H 2 0)*+ (IV) + c- -> An(H 2 0)^+ (III) . (1) 

Potentials are obtained from the free energies of the following three full reactions: 

An0 2 (H 2 0) 2+ (VT) + -H a (s) + H 2 -> An0 2 (H 2 0)+(V) + H 3 0+ 
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Experiment (A) Experiment (A) 



FIG. 2. Selected calculated mean internuclear distances compared to experimental values for (a) An=0 actinyl bonds and (b) 
An-0H2 bonds between the actinide atom and water molecules in the first solvation sphere. 




FIG. 3. Reduction potentials (volts) as calculated by DFT and corrected for the spin-orbit interaction compared to measured 
values as detailed in Ref. [45} see however Ref. |46|for a different compilation of measured redox potentials with slightly different 
values. PBE and B3LYP DFT results, corrected for spin-orbit and multiplet interactions, are from Table 10 of Ref. 24. 
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TABLE I. Calculated and selected experimental internuclear distances (A). The calculated lengths are means over all bonds of 
the given type. 



V 



VI 



u 



Np 



Pu 



Am 



U 



Np 



Pu 



Am 



An=O c 



Calc. 



Exp. 



1.84 



1.81 
1.8l[ a 

1.8. 



iff 



1.80 



1.79 
1.7C 
1.7S 
1.70 



1.76 1.74 
1.7^] 1.7^ 



1.74 



An-OH 2 



Calc. 
Exp. 



2.61 



2.60 



2.47 
2.5C 



2.62 
~2A^ 



2.65 



2.49 



2.47 



2.47 



2.4l[[ 

2.421* 



2.420 2.4C0 



2.50 



III 



IV 







U 


Np 




Pu 


Am 


U 




Np 


Pu 


An-OH 2 


Calc. 


2.53 


2.52 


2.50 


2.50 


2.40 


2.38 


2.38 


Exp. 


2.6lQ 


2.52 
2.51 


g 

c 


2.490 


2.4^] 


2.41 
2.4C 


C 

c 
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a Reference 1361 
b Reference [371 
c Reference 1381 

d References 391 and l40l as cited in Reference 1171 
c Reference 1411 
f Reference [42] 
s Reference 1431 
h Reference 1441 



An0 2 (H 2 0)+(V) + 3H 3 0+ + -H 2 (<?) -> An(H 2 0)^+(IV) + 2H 2 
An(H 2 0)^+(IV) + *H 2 (<?) + H 2 -> An(H 2 0)^+(III) + H 3 0+ 



(2) 



using the zero-potential reference half-reaction of the standard hydrogen electrode (SHE) H + + e~ — > |H 2 (g). We 
note that the free energy of the hydronium reaction H30 + — > H + + H 2 as calculated within DFT is —5.33 eV. 

Because Am(IV) is unstable, in the case of americium the V/III redox potential is calculated in lieu the V/IV and 
IV/III redox couples. In this case, the reduction half reaction is: 



Am0 2 (H 2 0)g"(V) + 4H 3 0+ + 2e" -> Am(H 2 0)g + (III) + 3H 2 
corresponding to the full reaction 

Am0 2 (H 2 0)+(V) + 2H 3 0+ + E 2 (g) -> Am(H 2 0)^+(III) + H 2 



(3) 



(4) 



Fig. [3] compares the calculated potentials to experiment; the ca lcula ted potentials include the contribution from 
the one-particle spin-orbit interaction as described below in Sec. Ill The VI/V potentials for U, Np, and Pu are 
qualitatively consistent with previous work despite differences in the DFT packages. We find redox potentials, in- 
cluding the spin-orbit interaction, of respectively —0.44, 0.31 and 0.94 volts. These compare to values 2.37, 4.00 and 
3.28 volts found by Hay, Martin, and SchreckenbacrP^ who used Gaussian 98, relativistic core potentials, an d the 
hybrid B3LYP functional, and included multiplet interactions corrections. However, Shamov and SchreckenbacrP^BS 
following a similar procedure but with a smaller (n < 4 frozen core obtained —0.10, 1.72 and 1.29 volts in better 
agreement with experiment (0.160, 1.236 and 0.966 voltsj^, highlighting the importance of treating all orbitals with 
n > 4 dynamically. Use of the Priroda PBE functional yielded comparable redox potentials of —0.51, 0.87 and 0.43 
voltsP3 When the multiplet interaction correction is removed these VI/V potentials they show the same trend as our 
calculated DFT potentials , increasing monotonically from U to Np to Pu. 

As mentioned above we only study the case of first solvation spheres of An(III) and An(IV) with 8 water molecules. 
Ref. |48]reports alternative molecular geometries for plutonium ions, with Pu(III) surrounded instead by 9 coordinating 
water molecules (see also Refs. [TBI and l2"Tj) . A DFT calculation of the free energy of Pu(III) with nine water molecules 
finds it to be lower by 0.09 eV, a small change in comparison to the other corrections that we consider here. 
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FIG. 4. KS molecular orbital energy levels and their atomic fragment composition for the Pu(IV) (top) and Pu(VI) (bottom) 
complexes, highlighting the importance of the 5f and axial (actinyl) 2p orbitals near the HOMO/LUMO boundary. The relative 
size of the contributions from the Lowdin atomic orbitals is indicated by the thickness of the lines. (Hydrogen orbitals are 
not shown for Pu(VI) and O includes contributions from both the actinyl ligands and the water molecules.) As indicated the 
atomic fragments have equal numbers of up and down electrons, but the complex has greater numbers of up electrons, as the 
spin-unrestricted DFT calculations satisfy Hund's rule. 



III. INDEPENDENT PARTICLE MODEL 



We turn next to the construction of the independent-particle model of the frontier orbitals that includes the spin- 
orbit interaction and serves as the starting point for the many-body Anderson impurity model: 

Hindep = Hks + H so . (5) 

We make the usual assumption that KS fermions may be regarded as physical electrons, when in fact the connection 
between the two is subtle: The DFT ground state is a Slater determinant of KS fermions whereas the electron 
wavefunction is not generally described by a single determinant. Conveniently ADF expresses the KS eigenstates 
as linear combinations of localized orbitals orthonormalized by the Lowdin procedure^. These Slater-type orbitals, 



G 



expressed in terms of Cartesian spherical harmonics, form the basis that we work in. In this way the KS orbitals are 
projected onto the Hilbert subspace that retains only the actinide 5f Lowdin orbitals and, in the case of the actinyls, 
2p Lowdin orbitals on the two oxygen ligands. As Fig. [4] shows, these are the atomic orbitals that are the most 
important constituents of the KS orbitals near the HOMO/LUMO boundary. For the actinyls, 14 + 6 + 6 = 26 Lowdin 
orbitals are thus retained; for oxidation states III and IV only the 14 5f orbitals are required. Although it would be 
desirable to also retain the An 6d Lowdin orbitals which are close in energy to the An 5f orbitals, and also the Lowdin 
orbitals in the first solvation sphere of water molecules, the resulting many-body Hilbert space would be too large for 
exact diagonalizations to be carried out. 

In the case of the actinyls the resulting effective Hamiltonian of the frontier orbitals may be written: 

H ks = 4j*fLfj<y + £ P ab*PLPb<i + ^{ticurfLPacr + H.C.) . (6) 
ija aba iaa 

Operator f! creates an electron in the 5f (2p) Lowdin orbital with spin a and spatial state i (a) in the Cartesian 

spherical harmonic basis. In Eq. J^ef.^ and e p abi7 are matrix elements for, respectively, the An 5f and O 2p Lowdin 
orbitals, and ti aa are the hopping amplitudes between the 5f and 2p orbitals. Oxidation states III are IV are modeled 
by the first term in Eq. [6] alone. Parameters e and t are obtained by calculating the matrix elements of the KS 
Hamiltonian, which is diagonal in the basis of KS orbitals, in the basis of the 5f (and in the case of actinyls, 2p) 
Lowdin orbitals. The matrix elements may then be grouped into the amplitudes that appear in Eq. [6] 

The truncation of the full Hilbert space to the subspace consisting only of An 5f and (actinyl) O 2p Lowdin orbitals 
introduces error into the calculations. Fig. [5] shows how the single-particle energies of the projected independent- 
particle model match those of the full unprojected KS levels for the cases of U(III) and U(V). In Fig. § the total 
electronic occupation of the An 5f and (actinyl) O 2p Lowdin orbitals as calculated by projection of the occupied KS 
molecular orbitals onto the Lowdin atomic orbitals is compared to the corresponding occupancies found by Mulliken 
analysis. The occupancy of oxidation states III and IV agree to within 2% of that found by Mulliken analysis; and to 
within 4% for oxidation states V and VI. 

As it stands Eq. [5] has no spin-flip processes, reflecting the limitation of DFT which is formulated in terms of 
separate densities of spin-up and down electrons. The one-electron spin-orbit interaction for the 5f Lowdin orbitals, 
H so = ■ s + Co, is therefore added to the independent-particle model. Energies C and Co are obtained from ADF 
calculations on isolated (gas-phase) actinide ions. The spin-orbit interaction splits the energies of the j = 5/2 and 
7/2 states by Ae = e 7/2 - e 5/2 = |C = 0-67 eV (U), 0.78 eV (Np), 0.87 eV (Pu) and 1.06 eV (Am). The interaction 
also shifts the energy of all 5f Lowdin orbitals upwards in energy by Co = 0.18 eV (U), 0.23 eV (Np), 0.26 eV (Pu) 
and 0.29 eV (Am). 



IV. MANY-BODY MODEL 



The Coulomb repulsion between 5f electrons is taken into account by adding the corresponding two-body interaction 
term to the Hamiltonian. The two-body interaction 

H ee = — ^ f Im 1 ,m 2 ,m 3 ,m4 /mj Ji /m 2 (J2 ^ m 3 2 CO 

Tn 1 +m 2 = TTl 3+ m 4 

is conveniently represented in terms of the Coulomb matrix elements 

I mi ,m 2 ,m 3 , m4 = Y F L c {L \mi,m i )c {L \m 3 ,m 2 ) (8) 

£=0,2,4,6 

where L is the tot al orb ital angular momentum of t wo 5f electrons, are the Gaunt coefficients^^ and F L are 



the Slater integrals! 52 ^ As discussed below in Sec. VI the F 2 , F 4 , and F 6 integrals parameterize the energetics of 
rearrangements of the electrons in the 5f shell and hence may be determined from spectroscopic data. F°, however, 
is sensitive only to the total number of electrons in the 5f shell and is expected to be highly screened. We treat it is 
the one adjustable parameter in the hybrid calculation. Like -ffmdep, the electrostatic repulsion terms are transformed 
into the Cartesian spherical harmonic basis for the purpose of numerical calculations. We performed numerical tests 
to verify that the spectrum of Eq. [7] reproduces published results for isolated actinide atoms^. 

As the interaction is already partially included at the DFT level, care must be taken to avoid double counting itP^ 
This we do by subtracting its contribution at the HF level of approximation, making an assumption, however, that 
DFT with the PBE functional is close to HF in its treatment of the interaction. Thus 

H = H in dep + H ee — H ee (9) 
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FIG. 5. Comparison of the KS molecular orbital energies against those of the independent-particle model that describes the 
projected 5f and 5f-2p subspaces, for the cases of U(III) and U(V). Shading of the KS orbitals indicate the relative size of their 
contribution to the 5f or 5f-2p Lowdin subspace, with the lightest lines representing orbitals that make little or no contribution. 

where the over line denotes the HF factorization of H ee . The one-body HF subtraction, H ee , that models the part of 
the Coulomb interaction already included in the DFT calculation is given by 



H e e „ ^ y i,Jmm' ^ram' (j) fmaf'n 



(10) 



m.m ,<7 



where the direct or Hartree interaction J is given by 



'mm' — / J -m,n,n' ,m 
n.n'a 



(fno-fn'v) 



(ii) 



and the exchange or Fock interaction K is 



K , -V/ 



m.n.m' ,n' 



(fnerfn'a} 



(12) 



The expectation values appearing in Eqs. 



11 



and 



12 



are calculated from DFT. By construction, then, (H ee — H ee ) = 



providing a valuable check on the numerical calculations. 
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FIG. 6. Total electronic occupation of the An 5f and (in the case of actinyls) O 2p Lowdin orbitals as calculated by projection 
of the occupied KS molecular orbitals onto the Lowdin atomic orbitals, compared to that found by Mulliken analysis. 



V. FRACTIONAL OCCUPANCY 

The combined occupancy of the 5f and 2p Lowdin orbitals is not an integer (see Fig. [6]). To handle this fractional 
occupancy within the many-body model of the reduced 5f-2p subspace, we calculate the ground state energies of the 
many-body model at the two integer occupancies that bracket the fractional value, and then compute a weighted 
average of the energies. The effectiveness of this algorithm may be illustrated with a simple model consisting of a 
single f-orbital and a single c-orbital, where the c-orbital is a model for orbitals not included in the restricted 5f-2p 
subspace. The model is parameterized by on-site f-orbital energy e, a hopping amplitude between the f- and c-orbitals 
t, and the Coulomb repulsion U between two electrons in the f-orbital: 

H = eflf a -t (/t c a + H.c.) + U flflhfr (13) 

where the sum over the repeated a spin index is implied. The model is easily diagonalized. For instance in the 
2-particle, spin-singlet, subspace spanned by the 3 basis vectors: 

|l)=4c{|0) 



|2>^(44-ct/t)|0) 

|3> = /}/>) (14) 

H takes the form of a 3 x 3 matrix and, when diagonalized, the resulting exact ground state energy E may be 
compared against approximations. 
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TABLE II. Comparison of ground state energies for the two-orbital model with t = 1 and e = —3. Here (n) is the occupancy 
of the f-orbital in the HF approximation; Eo is the HF energy; Eo is the improved estimate of the ground state energy; and Eo 
is the exact ground state energy. 



u 


(n) 


E 


E 


Eo 


4 


1.579 


-3.874 


-4.051 


-4.323 


6 


1.407 


-3.079 


-3.607 


-4 


8 


1.246 


-2.571 


-3.708 


-3.860 



The HF approximation to the Hubbard interaction, H ee = Uf^fJ^fifo = Un^n^, is given by setting I = U in Eq 
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and using the fact that (n-j-) = (nj.) = \{n) by spin-rotational invariance of the spin-singlet ground state where 
n = n t + n± = f\fo + /j/j,. The result is: 

H ee = j(n)n (15) 

and it replaces the two-body interaction with one-body term that renormalizes the one-body f-orbital energy: 

e^e r ({n)) = e+j{n) . (16) 

Self-consistency is then attained by adjusting (n) so that the f-orbital occupancy as calculated in the ground state of 
the independent-particle Hamiltonian h: 

h = e r ((n)) ft U-t&c. + H.c) (17) 
equals (n) . The resulting HF equation is 

2 (18) 



l , ( ^((»)) + V^({"»+4« ; 



and it can be solved by iteration. The HF approximation to the ground-state energy of the two-electron system is 
then given by: 

E = e r ((n))- v/e2«n»+4i2 (19) 

which is simply the energy from filling the lowest eigenstate of the renormalized independent-particle model with both 
a spin-up and a spin-down electron. 

An improved approximation of the ground state energy can be obtained by carrying out an exact diagonalization 
in the reduced subspace consisting of only the f-orbital. Double-counting the interaction is avoided by subtracting 
the HF contribution to the Coulomb energy, H ee from the exact two-body Hubbard term Un^n^ = (U/2)(n 2 — n). 
The improved estimate of the ground state energy is given, for fixed integer occupancy n = 0, 1, 2 of the f-orbital, by 
Eq. [9] which in this simplified context reads: 

(n) =E + ^ (n 2 -n)-j(n)n. (20) 

Finally a weighted average of Eo(n) based on the HF occupancy (n) = 1 + x yields, for x > 0: 

E = (1 - x) E (l) + x E Q (2) (21) 

and as Table |ll| shows, for the case of t = 1 and e = —3 there is a substantial improvement over the HF approximation. 
The ground state energy decreases because the two electrons are now correlated and able to avoid each other. 



VI. RESULTS 

The Slater integrals F 2 , F 4 , and F 6 parameterize changes in electrostatic energy due to rearrangements of the 
electrons in the 5f shell, and as a consequence are insensitive to the chemical environment surrounding the actinide. 
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TABLE III. Slater integrals (eV) used in the Anderson model. 



Ion Configuration F 2 F 4 F fi 

U 4+ 5f 2 5.746 3.693 2.201 

Np 4+ 5f 3 6.249 4.016 2.394 

Pu 4+ 51 4 6.778 4.356 2.597 

Am 4+ 5f 5 7.867 5.056 3.014 



We use the values displayed in Table III These are based upon spectroscopic data^ that is expressed in terms of 



Racah parameters E , E 2 , and E 3 . The Racah parameters are linearly related to the Slater integrals by the following 
equations: 

530^ 
9 

396^3 

^ • (22) 



F 2 = 




6895 ~ { 
9 E ~ + - 


F 4 = 




12870 o 

-^-E 2 + 


F 6 = 


9 


3904 - 
f — E 2 ~ 



The values for F 4 and F 6 listed in Table III are somewhat lower than those obtained from Ref. [56] (a range of 
values may be found in the literature, see for instance Ref. l54|) but we have checked that the corrections to the redox 
potentials are insensitive to these differences. By contrast F° parameterizes the part of the Coulomb energy that is 
sensitive to the total number of electrons occupying the 5f shell and is thus important in charge-transfer reac tions 



As expected the corrections to the redox potentials vary with F° (see below). As discussed above in Sec. IV F° is 
expected to be highly screened, but the degree of screening is difficult to predict reliably from first-principles. We 
therefore treat it as the one adjustable parameter in our calculations. 

The many-body Hilbert space has a maximum dimension 26!/13! 2 = 10,400,600 for the case of 13 electrons 
populating the 5f and 2p Lowdin orbitals of an actinyl. Exact diagonalization of the many-body Hamiltonian is 
accomplished with the use of the sparse-matrix Davidson algorithm. As the HF approximation may be formulated 
as a variational problem over the subset of wavefunctions that are single Slater determinants, in the absence of the 
spin-orbit interaction diagonalization of H, Eq. [9J yields ground state energies that are less than the ground state 
energy of Hi n de P , Eq. tal The reduction in the energy is a consequence of the fact that correlations between the 5f 
electrons permit the electrons to avoid each other more effectively than when the interaction is described only at 
the mean-field level. The effect, which holds even in the presence of the spin-orbit interaction, is evident in Fig. [7] 
where it can be seen that the electronic ground state energy of the plutonium complexes decreases with increasing 
F°. Also as F° increases, electrons in the actinyls move from the 5f Lowdin orbitals to the 2p Lowdin orbitals of 
the oxygen ligands (the occupancy remains fixed for Pu(III) and Pu(IV) because only the 5f orbitals are included in 
the many-body model). This charge transfer is also evident in the total spin in the 5f orbitals which decreases in the 
actinyls with increasing F°. These observables, like the ground state energy, are calculated from a weighted average 
of the two many-body ground states of integer occupancies that bracket the occupancy obtained from DFT. 

The many-body correction to the ground state energy as computed within DFT changes the free energies of the 
reactions Eqs. [2] and [4] and hence the redox potentials. Fig. [8] shows the change in the redox potentials as a function 
of F°. Potentials for F° — are the DFT values corrected by the spin-orbit interaction and the Slater integrals F 2 , 
F 4 , and F 6 . For the VI/V couples (shown in blue), the closest match with experiment is for F° = 2.1 eV (U), 1.3 eV 
(Np), 0.4 eV (Pu) and 3.7 eV (Am). For the IV /III redox couples (red), the best match occurs for F° = 1.4 eV (Np) 
and 0.1 eV (Pu) eV, but the correction to the electronic energies worsens the match to experiment in the case of U. 
However, in the case of the V/IV redox potentials (green), as well as the Am V/III potential (black), the pure DFT 
values are all low in comparison to experiment (see Fig. [3]), and the corrections only lower the potentials further. 

As it stands the calculation does not account for changes in the screening of the Coulomb interaction as the oxidation 
state changes. A reduction in the size of F° for the actinyls yields in most cases a better match with experiment, 
particularly in the case of the V/IV redox couple. In Fig. ^F° = 2.5 eV for U(III) and U(IV); 1.0 eV for Np(III) and 
(IV); 0.5 eV for Pu(III) and (IV); and 0.5 eV for Am(III). In each case F° is reduced by 0.5 eV for oxidation states V 
and VI. There is a striking improvement in the match with experiment with two exceptions: The calculated U IV/III 
and Am VI/V potentials move further away from the measured values. The downward shift in F° may reflect the role 
that the oxygen ligands play in screening repulsion between 5f electrons. A similar trend has been found in studies 
of isolated molecules and ions. For example in Ref. [S3] a range of values for F° — 2.3 to 3.3 eV are reported for an 
isolated U 4+ ion, decreasing to F° — 1.6 eV in the case of the UPt3 molecule. 
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FIG. 7. The correction to the DFT ground state energy, the occupancy of the 5f Lowdin orbitals, the total spin (S 2 ) in the 5f 
orbitals, the total orbital angular momentum (L 2 ) in the 5f orbitals, and the total angular momentum (J 2 ) in the 5f orbitals 
as F° is varied, holding F 2 , F 4 , and F 6 and the spin-orbit interaction fixed. The four oxidation states of plutonium are shown. 



Fig. 10 presents probability distributions of different electron occupancie d 57 * 58 ! in the actinyl 5f Lowdin orbitals, 



using the same values of F° adopted in Fig. [9] The distributions are particularly broad for neptunium and plutonium 
reflecting the increasing number of 5f electrons as one moves along the row of early actinides competing against an 
increasing tendency to localize. 
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FIG. 8. Calculated redox potentials of the actinide redox couples An(IV)/An(III) (red), An(V)/An(IV) (green), An(VI)/An(V) 
(blue) and Am(V)/ Am(III) (black) as a function of F° . Experimentally determined potentials are indicated by the solid 
horizontal lines. 



VII. CONCLUSION 



The hybrid DFT / many-body approach taken in this paper shows some promise for the theoretical modeling of the 
difficult but important problem of redox chemistry involving the early actinide elements. By incorporating the physics 
of strong correlations between electrons in the frontier orbitals we are able to correct the electronic contribution to 
the free energy as computed in DFT, and thereby bring the calculated redox potentials into closer agreement with 
measured values. The calculations require one adjustable parameter, F° , for each actinide species, in addition to a 0.5 
eV downshift in F° for the actinyls, yet has predictive power as it yields potentials for 3 redox couples An(VI)/An(V), 
An(V)/An(IV), and An(IV)/An(III). 

The match with experiment is improved in 6 of the 11 redox reactions that we study; agreement remains good in the 
case of 2 reactions (Pu VI/V and Np IV/III), little changed but poor for Np V/IV, and worsens for 2 others (U IV/III 
and Am VI/V. The calculated potentials are certainly not of chemical accuracy, but they do evidence significant 
trends. In the case of plutonium, for instance, the calculated potentials approach the measured near-degeneracy of 
the 3 redox potentials, and hence go some distance towards explaining the propensity of plutonium species in solution 
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FIG. 9. Many-body corrected redox potentials (filled symbols) compared to measured values. Open symbols are the potentials 
obtained from DFT including only the spin-orbit correction. F° = 2.5 eV for U(III) and U(IV); 1.0 eV for Np(III) and (IV); 
0.5 eV for Pu(III) and (IV); and 0.5 eV for Am(III). In each case F° is reduced by 0.5 eV for oxidation states V and VI to 
account for increased screening of the Coulomb interaction in the actinyls. 



to easily disproportionate and co-exist in several different oxidation statesP In the language of Hubbard models, 
disproportionation may be viewed as a consequence of an effective negative-U interaction for the overall complex; see 
for instance Refs. ">!) (TT1 In our calculations its origin may be traced in part to the strong correlations between the 
5f electrons that permit the electrons to avoid each other more efficiently than they can at the level of LDA/GGA 
or HF, lowering the electronic energy. When combined with all the other contributions to the free energy (solvation, 
vibrations, and translations) the redox potentials become degenerate and an overall effective negative-U interaction 
emerges. 

The calculations can be improved or extended in several ways. At the DFT level, different geometries with varying 
numbers of water molecules in the solvation spheres can be investigated. It may be desirable to treat a second sphere 
of solvation quantum mechanically rather than with continuum dielectri c mo dels. Of geochemical interest is actinidc 
complexation with carbonate and silicate substrates, and with colloids P^2I anc j these could be investigated by the 
hybrid approach. The projection onto Lowdin orbitals could be replaced with projective orthogonalizatiorp^ to better 
minimize the admixture of neglected orbitals. It may be possible to include some additional orbitals in the many-body 
model, if not by exact diagonalization then possibly by methods such as the density-matrix renormalization-group 
(DMRG). It would also be interesting to investigate the problem of double-counting the interaction in the hybrid 
approach by carrying out a pure HF calculation, working directly with the physical electrons rather than with KS 
fermions. However it is expected that the pure HF calculation will not by itself be a sufficiently accurate foundation 
for free energy calculations. Ultimately it would be desirable to replace hybrid DFT / Anderson impurity model 
approach developed here with a unified first principles method that can simultaneously and accurately describe the 
physics and chemistry of relativity, solvation, and strong electronic correlations. 
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